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Abstract 

A remarkable property of Hermitian ensembles is their universal behavior, that is, once properly 
rescaled the eigenvalue statistics does not depend on particularities of the ensemble. Recently, 
normal matrix ensembles have attracted increasing attention, however, questions on universality 
for these ensembles still remain under debate. We analyze the universality properties of random 
normal ensembles. We show that the concept of universality used for Hermitian ensembles cannot 
be directly extrapolated to normal ensembles. Moreover, we show that the eigenvalue statistics 
of random normal matrices with radially symmetric potential can be made universal under a 
conformal transformation. 
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One of the most celebrated results in random matrix theory is the universal behavior of the 
eigenvalues statistics of Hermitian ensembles. The universality states that the eigenvalues 
statistics considered on the proper scale (with unit the mean eigenvalue spacing) is always 
the same, independently on the particularities of the ensemble. The universality in random 
matrix theories plays a major role in many branches of sciences l[ [j, [3, 4]. 

The universality problem is addressed by using the so-called invariant ensemble model. 
The probability 

P(M)alM cxexp{-nTr[F(M)]}dM, (1) 

of finding a n x n matrix M within the ensemble, with alM the Riemann volume in the space 
of normal matrices, is invariant by unitary transformations. The corresponding eigenvalue 
density, in the limit n — > 00, depends on the particular form of V(M). 

The universality in the eigenvalue statistics can be heuristically understood in terms of 
the Dyson interpretation of the eigenvalues as a Coulomb gas For a hermitian matrix 
Eq. ([I]) can be written in terms of the joint probability of the eigenvalues 

n 

P„(Ai,--- ,A n )ocexp{ Yl M^-Ajl^-nJ^CAi)}. (2) 

l<iy^j<n i=l 

The joint distribution of eigenvalues has the same form as a 2 dimensional Coulomb gas 
restricted to one dimension. The eigenvalue interaction is given by a Coulomb term In |Aj — 
A j I — 1 under a potential V(X) that confines the eigenvalues in a bounded set of the real line. 

Since the eigenvalues are restricted to one dimension, if an eigenvalue is located between 
two eigenvalues the Coulomb potential does not allow it to switch its position with the 
surrounding eigenvalues. As the number of eigenvalues increase (the size of the matrix) the 
eigenvalues become closer and the effect of the Coulomb potential overtakes the potential 
V(X) (due to the logarithm behavior near zero). 

If we then rescale the mean distance between eigenvalues to the unity, heuristically this 
is equivalence to switch off the external potential. Thus no matter the external potential, 
after the scaling only the Co„,o,„b potential affects the eige„va.ue, Reeett.y, all this 
reasoning has been rigorously proven (see e.g. [5| and references therein). 

Over the last decade, normal matrix ensembles have attracted a great deal of attention, 



mainly due to their broad applications such as in quantum Hall effect 







quantum flows 



quantum field theory 8|, conformal hierarchies and conformal maps [9, 1Q|. Despite of 
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the increasing interest on normal ensembles, questions on their universality properties still 
remain open. 

The joint probability of the eigenvalues of a normal ensemble obeys an equation of the 
type of Eq. (j2J). The basic difference is the dimension of the eigenvalues Zj. For normal 
matrices the eigenvalues lie on the complex plane. The different eigenvalue dimensions 
between Hermitian and normal ensembles may lead to different universal behaviors. In the 
later case, as we add eigenvalues they can turn around and occupy a variety of positions in 
the complex plane to finally be set to minimize the electrostatic energy. It is possible that 
after this process the rescaling of the average distance still brings out the particularities of 
the potential. If so, what would be the procedure to assure the universal behavior in normal 
ensemble, if any? 

In this letter, we show that ensembles of normal matrices under radially symmetric po- 
tential are not universal in the Hermitian sense, that is, if one rescales the mean distance 
to unity the eigenvalue statistics still depends on the potential. We then show that the uni- 
versal properties can be retrieved by conformal transformations of the rescaled eigenvalue 
distribution. 

Recently, questions on universality of normal ensembles have been addressed in Ref. 11]. 



The machinery of statistical mechanics was used to argue that normal ensembles presents a 



universality of the second type 



1 2t ] . This means that the one-point Green's function and the 



connected two-point Green's functions may be related to each other in a universal way after 
some proper normalization and rescaling of the system. Yet, such proper rescaling remains 
unknown, and even more serious question remains: whether universality of the second type 
is the only kind of universality ruling normal ensembles 13]. 

Hermitian ensembles: Let us first revise the precise meaning of universality for Hermitian 
ensembles. The quantity playing a major role in the analysis of random ensembles is the 
integral kernel 

n 

K n (x, y) = e -SCVM+V(y)) £ ^ {x) ^ {y) (3) 

3=1 

where V(x) is the potential [see Eq. (j2J)] and {4>j}™ =1 denotes the set of polynomials 
up to order n — 1, orthogonal with respect to the weight exp{— nV(x)}. The important 
statistical quantities associated with the matrix ensemble such as the eigenvalue density 
jo(Ai) = J P n (Xx, • ■ ■ , \ n )d\ 2 • • • A n can be obtained by the K n ^]. For large n the relation 
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reads 

K n (X,X)=np(X)(l + o(l)), 

where o(l) converges to as n — > oo. 

To proceed the universality analysis, we must rescale the integral kernel. The scale is 
chosen such that the fraction of eigenvalues in an interval of length s close to a point A 
equals s, in other words, the average spacing between eigenvalues is unity. It can be shown 
that the correct scale is K n (X,X). Indeed, since the fraction of eigenvalues in the interval 
i C lis given by 

1 n 

f n (A) = -X)xA(Ai), (4) 

i=i 

it is possible to show that, for large n, it holds (nf n [X, X + s/K n (X, A)]) = s[l + o(l)] [14j . 
Once we have the proper scale K n (X,X), we proceed the analysis by rescaling the integral 
kernel 



The astonishing result [14] is then 

~ sinvr(x-?/) 
hm K n [x, y) = — (6) 

n-s-oo 71{X — y) 

exists pointwise for Freud-type or real analytic potentials V. Since it does not depend on 
V the Hermitian ensembles in this sense are universal. The natural questions is whether 
normal ensembles also display universal eigenvalue statistics. 

Radially symmetric normal ensembles: We shall adapt these questions on universality 
to the family of normal ensembles. We shall focus on the case of normal ensembles with a 
radial potential 

V a (z) = \z\ a . (7) 

we shall consider the potential for a > 2. The statistical quantities of this ensemble can 
be characterized b y u sing the theory of logarithmic potentials and orthogonal polynomials 



developed in Ref. [15|]. The eigenvalue density associated with the potential V a (z) is given 
by p(z) = a 2 \z\ a ~ 2 /47r, the support supp(p) of the eigenvalue density can be computed by 
the condition |,2|V^(|z|) < 2, which yields 

supp(p) = {zeC:\z\< (2/a) 1/a }. 



The integral kernel associated with V a is then given by K%(z,w) = exp{— |[V^(z) + 



(w), where (j)j(w) are orthogonal polynomials with respect to the 



weights exp{nV a (z)} . Since, V a (z) is radially symmetric the orthogonal polynomials <t>j{w) 
are monomials, using the explicit relations for <pj(z) we obtain 

_»L,,|a n 2j. 
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2tczw 



n 



(8) 



For n large enough K®(z,w) can be made small everywhere, except in a vicinity where 
K®(z, w) attains its maximum, which occurs when z = w. In Fig. ([!]), for z = 0.3 + iOA we 



depict n 1 



KZ(z,w) 



as a function of w = x + iy. One can clearly see that as n increases 



the peak is sharped at w 
of z e supp(p). 



z. n 



K%(z,w) 



presents the same behavior for different values 




FIG. 1: Profile of the integral kernel n _1 K®(z,w) for fixed z = 0.3 + i0.4 as a function of 
w = x + iy. We show distinct values of a and n. One can see that the kernel presents a peak that 
is sharped with the increase of n. This behavior persists for all a > 2. The distinct values of a are 
a = 4 in (A,B) and n = 50, 200 and a = 8 for n = 50, 200 in (C,D). 



Asymptotic expansion for the integral kernel: We wish to obtain an asymptotic expression 
for the integral kernel K®. We are concerned with the behavior of the integral kernel for 
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all w in a vicinity of z, only the values of w close z gives a relevant contribution. We shall 
investigate the behavior of K% in the region 

A a := |z, w G supp(p) : |argz — argu>| < — j, 

where argz stands for the argument of the complex number z. We shall explain this partic- 
ular choice latter on. For now it is important that A a contains the peak of K"(z, to). 

In the vicinity A a we can approximate Eq. (jHJ) by an integral and then solve it by means 
of the steepest descent method. Proceeding in this way, we obtain an asymptotic expansion 
for the integral kernel. The technical details of this analysis shall be given elsewhere. The 
asymptotic expansion reads 

n- 1 k%(z,w)(l + o(l)) if(z,w)£A a 
n- 1 K%(z,w) = \ '\ >) (9) 

otherwise 

i. 

with the leading part of the asymptotic expansion being 

-K%(z,w) = — (zW)* -1 exp|n[(zW)T - — — ]|. 

Now it becomes clear the behavior seen in Fig. (0Q). Since 2Ke(zw) a ^ 2 < \z\ a + \ w\ a , being 
equal only if z — w, with the increase of n, K® converges to zero exponentially fast whenever 
z 7^ w. 

The o(l) behavior of the asymptotic expansion can be seen by introducing the quantity 



= max 

z,w£supp(p) 



n- x K a n { Z ,w)-n x K a n {z,w) 



which gives the largest value of the difference between the two kernels in A a . In Fig. ([2]) we 
depict R® as a function of n for different values of a. One can clearly see the behavior o(l) 
ofiC 

Outside the region A a the quantity n K® may differ from n~ l K". \K"(z,w)\ has a 
term exp{n(zw) a ^ 2 } , whose absolute value is a periodic function of period of the argument 
(zw). To see this we write again z = re ld and w = se % ^, which yields 

\exp{n(zw) a/2 }\ = exp{n(rs) Q/2 cos[a(#-0/2]}. 

As a consequence other peaks may occur in an angular distance 4n/a far apart from the 
main peak localized in A a . 
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FIG. 2: We depict the maximum difference versus n. We plot i?" f° r three distinct values of 
a, namely a = 6 (circles), a = 8 (squares) and a = 11 (diamonds). One can clearly see the o(l) 
behavior of 

Rescaled integral kernel: As for the Hermitian ensembles, we shall rescale the average 
spacing between eigenvalues. Denoting by D s the disk of radius s in the neighborhood of 
the origin, we must find a scale g(s) such that the fraction of the eigenvalues in D s equals 
the area of the disk. We must have (nf n (D g ^)) = ns 2 [l + o(l)], with /„ the indicator 
function as defined before. Explicitly, for large n through the equilibrium density pire ), 
using z = re ld we have 



' Jo Jo 

solving the integral it yields ((g(s)) a ) (1 + o(l)) = 7rs 2 [l + o(1)]. The scaling function has 
the form 



Considering only terms of superior order in n the rescaled integral kernel is given by 



Note that the difference between K a (z,w) and the rescaled kernel for the Hermitian 





(10) 



K a (z, w) = g' (z) K« (g (z) , g (w)) g> (w). 



(11) 
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ensemble in Eq. (E]) is due to the Jacobian of the coordinate change of w h g{w) in the 
complex plane. After some manipulations the rescaled kernel reads 

K a (z,w) = -e a V 2 2 ). (12) 
a 

Hence, the rescaled integral kernel, unlike the rescaled kernel for Hermitian ensembles, 
still carries peculiarities of the potential. Therefore, the concept of universality used for 
Hermitian ensembles cannot be directly extrapolated to normal ensembles. The question 
then turns to whether there is a scheme to obtain a universal integral kernel and consequently 
a universal eigenvalue statistics. 

Conformal Universality: In the following, we shall show that there is a conformal map <p 
such that 

K a k 

where K does not depend on the potential V a and, in this sense, is a universal kernel. The 
universal kernel 

1 - H 2 M 2 

K{z,w) = -e zw -— -— . (13) 

7T 

The conformal map ip : £ a (z) — > C reads 

(p{w) = \pnw^ . 

where T, a (z) = {w G C : |arg2; — argw| < 27c/a}. The explicit relation between of integral 
kernels reads 

K n (z, w) = ip'{z)K{ip{z), (^H)TH- (14) 
Since ip is a conformal transformation, we have 

K(z,w) ^ K:(z,w) 

where <p~ l (w) = (w/n) 2 ^ a maps C into S Q (z). Thus, we conclude that the family of poten- 
tials V a (z) = \z\ a is conformally universal. 

It is remarkable the fact that (p~ l is intimately connected to unitary rescaling g. Intro- 
ducing u(w) = ^a/2TTW, we may write 

g{w) = (cp^ 1 ou^iw), 
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where u(w), acting as a rescaling of the complex plane, may stretch or shrink the distances 
depending on a. Hence, the unitary rescaling g is associated with the universal behavior of 
the ensemble after a proper rescaling of the complex plane. 

Our results have also shown that for normal ensembles there is a distinction between K a 
and the universal kernel K(z,w). Remarkably, the relation between K and K a is realized 
by a rescaling of the complex plane through the conformal map u(z) = y/l^z. The explicit 
relation reads 

K a (z, w) = u'(z)K(u(z),u(w))u'(w). (15) 

The difference between K and K a is not observed in the Hermitian case. In the normal 
case the difference can be heuristically explained by the difference in the dimension of the 
eigenvalues. For the Hermitian case the influence of the Coulomb repulsion seems to be 
much more relevant than for the normal ensemble. Moreover, the equilibrium disposition 
of the eigenvalues varies with the parameter a due to the many possible stable crystalline 
structures. 

In summary, we have shown that the concept of universality used for Hermitian ensembles 
cannot be direct extrapolated to normal ensembles. We have considered a model for random 
normal ensemble, where the potential is given by a radially symmetric function depending 
only on one parameter. We have shown that the rescaling of the average distance between 
the eigenvalues to the unity, unlike for Hermitian ensembles, does not lead to a universal 
eigenvalue statistics. To overcome this difficult, we have put forward a new concept of 
universality for normal ensembles, which we called conformal universality. We have shown 
that the rescaled integral kernel K"(z,w) can be obtained by conformal transformations 
from a universal kernel K(z,w). Our result shows that the intricate relationship between 
conformal geometry, revealed in the recent works [9J, HOj , might play an even more important 
role than previously thought. It is wishful to extend the procedure adopted here to general 
potentials. 
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detailed reading of the manuscript. We would like to acknowledge the financial of the 
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